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Abstract 

We study the process of equilibration between two non-extensive subsystems in the framework 
of a particular non-extensive Boltzmann equation. We have found that even subsystems with 
different non-extensive properties achieve a common equilibrium distribution. 
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1. Introduction 

Power-law tailed spectra occur in a wide range of physical phenomena, re-association 
in folded proteins [1], fluxes of cosmic rays [2], turbulence [3], finance and economics [4], 
electron-positron annihilation [5] , motion of Hydra cells [6] , epilepsy [7] , linguistics [8] , 
nuclear physics [9], astrophysics [10], field theories and cosmology [11], scientific citations 
[12], distributions of individual succes of musicians [13], urban agglomerations [14], in- 
ternet phenomena [15], phase transformations [16], in algorithms for global optimization 
(e.g., the generalized simulated annealing) and related computational methods [17] and 
information theory [18]. 

We are particularly interested in heavy ion collisions where cut power-law distributions 
describe transverse momentum (j>t) spectra at low and intermediate values [19]. The con- 
ventional approach [20] to the distribution function fits a Boltzmann- Gibbs exponential 
characterized by a single temperature and the deviation from the exponential at the tail 
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is attained to non-equilibrium effects. A transverse expansion distorts the locally thermal 
spectra; at low px by suppression, at high pt by a blue shift factor in the temperature. 

Some further attempts have been made to interpret power-law tailed spectra as equi- 
librium phenomena for the whole pr range: Non-extensive thcrmodynamicses predict 
such distributions [9,21]. It has been shown, that dynamicses among (quasi-) particles 
can be constructed, leading to power-law tailed stationary distributions [22]. This parton 
cascade model, referred to as non-extensive Boltzmann equation (NEBE), generates sta- 
tionary states subject to non-extensive kinetic energy addition rules. We have shown that 
an infinite class of equilibrium distributions can be established depending on the way we 
generalize the Boltzmann equation [23]. In this paper we point out that cut power-law 
distributions play a special role among all possible functions: The rule applied in the 
non-extensive Boltzmann equation leading to such distributions is the leading term in a 
low-energy expansion of a general associative rule. 

A crucial question towards non-extensive thermodynamics has been discussed in [24] , 
namely how do two non-extensive thermodynamical systems thermalize, if at all. In this 
article we aim to answer this question by performing momentum-space simulations within 
the NEBE model. 

2. Non-Extensive Boltzmann Equation (NEBE) 

In the non-extensive extension of the Boltzmann equation we keep the original factor- 
izing form for the two-body distributions, 

jrji = f W1234 C/3/4 - /1/2) , (1) 

01 J234 

but in the transition rate 

^1234 = Mf^ • S (( Pl + ft) - (p 3 + p 4 )) • 5 (h (E u E 2 ) - h (E 3 ,E 4 )) , (2) 

the two-body energy composition is generalized to a rule h (Ei,E 2 ), which is not neces- 
sarily the simple addition. There can be many physical sources for this deformation of the 
kinetic energy addition rule; the most general being a pair interaction due to a potential 
whose value differs before and after the two body collision. This way corrections occur, 
which (e. g. by using a virial theorem) may be expressible in ter ms of the individual 
kinetic energies of the colliding subsystems (particles). Here E{ = y/pf + mf are the free 
kinetic energies of relativistic particles with mass mi. In a heavy- ion collision they are 
regarded as the asymptotic energies detected after desintegration of the system. 

The energy composition rule h(Ei, E 2 ) contains contributions stemming from the pair 
interaction. It is not trivial whether these can always be divided to one-particle contribu- 
tions, supporting a quasi-particle picture. In the followings we demonstrate that under 
quite general assumptions about the function h(x, y) the division of the total energy 
among free particles can be done. 

We assume that the generalized energy sum is associative, 

h(h(x,y),z) = h(x,h(y,z)). (3) 

Then due to a mathematical theorem [25] a strict monotonic function X(h) maps the 
energy composition rule to additivity 

X(h)=X(x)+X(y). (4) 
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This solution (4) of the functional equation (3) is unique up to a constant multiplicative 
factor. In this case the stationary solution of the non-extensive Boltzmann equation is 
given by 

/(P) = ^- X(£)/T , (5) 

and X{E) is regarded as the energy of a quasi-particle. 

One finds several examples for the energy composition rule h(x,y). Some of them 
lead to thcrmodynamicses where the entropy formula or the distribution function is 
familiar: e.g. h(x,y) = axy leads to the Renyi-entropy Sr — In J f 9 , with s = 

j-— — , q=l- aT, and X(E) = \ \n(aE). The distribution function is given by f(E) = 

±{aE)-V( aT \ 

The energy composition rule 

h(x, y) = x + y + axy (6) 

has the solution 

X(E) = -\n{l + aE), (7) 
a 

and leads to the cut power-law stationary distribution 

f(E) = 1(1 + aE)- 1 / aT . (8) 

(cf. for q = 1 — aT the q— Tsallis distribution emerges [23]) 

This form is the next to leading order expansion of a general associative rule for low 
energy: By Taylor expanding eq. (4) for low x, y and h values and requiring h(x, 0) = x, 
h(Q, y) — y one arrives at 

W x *"(0) 

h(x,y)=x + y-—-xy + .... (9) 

We conclude that the Tsallis-type composition rule with a = —X"(0)/X'(Q) is generic 
for leading order non-extensive effects at low energy. 



3. Non-extensive parton cascade 



In order to investigate the equilibration of non-extensive systems we start with two 
subsystems, equilibrated separately. In order to prepare these systems the non-extensive 
Boltzmann equation is solved numerically in a parton cascade simulation [22]. We use 
several different initial momentum distributions, and then make random binary collisions 
between randomly chosen pairs of particles. By doing so we apply the rules 

X(E l )+X(E 2 ) = X (E 3 ) + X (E 4 ) , (10) 



Pl+P 2 =P3+P4- ( n ) 

In each step of the simulation we select two particles to collide. Then we find the value for 
the new momentum of the first particle (p 3 ) satisfying the above constraints but otherwise 
random. Then applying eq. (11) we calculate the momentum of the second outgoing 
particle (p 4 ). In these particular simulations, this paper reports about, we use the free 
dispersion relation for massless particles (Ei(p i ) = We proceed with the next 
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collision repeating the above steps with a new randomly chosen pair of particles. A typical 
simulation includes 10 6 — 10 7 collisions among 10 5 — 10 6 particles. After 3 — 5 collisions 
per particle on the average, the one-particle distribution approaches its stationary form 
sufficiently. 

The following quantities are conserved during the simulation: 

N N N 

x{E tot ) = Y J xm, p = y.P" N = T, 1 - ( 12 ) 

i— 1 i— 1 i— 1 

We use eq. (6) for the energy composition rule, where a is called the non-extensivity 
parameter. Our model reconstructs the traditional Boltzmann-Gibbs thermodynamics in 
the limit of a = 0. More details can be found in [22]. 

We perform simulations on different systems with particle numbers Ni and -ZV 2 , total 
(quasi-)energies X(E\) and X{E2) and non-extensivity parameters a\ and a-i- We evolve 
these systems until they reach their stationary states. The unified system is taken as an 
initial state with N = N\ + N2 particles. 

Now we have three non-extensivity parameters 01, 0,2 and ai2, describing the three 
possible interaction types in the composed system. a\ is the parameter of the collisions 
among particles in the first subsystem, similarly a 2 corresponds to the second subsystem, 
and finally an is the non-extensivity parameter in a collision between a particle from 
the first subsystem and a particle from the second subsystem. We investigate cases with 
common and different a parameters. 



4. Results 




Fig. 1. a, Equilibration of two Boltzmann-Gibbs systems (a± =0,2= a\2 = 0) ; b, Equilibration of two 
Tsallis-type non-extensive systems (ai = (22 = a\2 = 2) 

Our results show that the subsystems equilibrate, in the final state of the composed 
system they have a common stationary distribution. 
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Fig. 2. a, Equilibration of two Boltzmann-Gibbs systems (ai = <22 = 012 = 0) ; b, Equilibration of two 
Tsallis-type non-extensive systems (ai = (22 = a\2 = 2) 




Fig. 3. a, Equilibration of a Boltzmann-Gibbs and a Tsallis-type non-extensive system (ai = 0, 02 = 2 
and ai2 = 0) ; 6, Equilibration of a Boltzmann-Gibbs and a Tsallis-type non-extensive system (ai = 0, 
02 = 2 and 012 = 2) 

We present examples with different initial conditions. In all of these simulations the 
particle numbers are the same for each subsystem, N\ — N2 — 250 000. The number of 
collisions in one simulation is N co u = 5 000 000, so N co u / (N1+N2) — 10 collisions happen 
per particle. Figure 1 shows the average energy and quasi-energy versus the number of 
collisions per particle. On each part of the figures 2 and 3 four curves show the initial 
(continuous line) and final (dashed line) energy distributions of subsystems 1 and 2. 
It is hard to distinguish the final states of the subsystems, because the corresponding 
distributions are very near to each other. Within numerical uncertainities they have a 
common stationary energy distribution. 

Figures la and 2a show the equilibration of two Boltzmann-Gibbs systems, where a\, 
<i2 and a\2 are equal to 0. This simulation was done for test purposes. As expected, in 
the final state the two subsystems have a common stationary energy distribution with a 
common temperature. The energy content per particle of the composed system became 
the arithmetic mean of the respective initial values. On figure la the curves for average 
X av — X(Ei)/N coincides with the curve for E av = J2i Ei/N, which comes from the 
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definition of the quasi-particle energy (cf. eq. (7) with a = 0). 

The equilibration of two Tsallis-type non-extensive systems with a\ = a 2 = an = 2 is 
shown on figures 16 and 26. The initial quasi-energy contents are different but at the end 
this difference disappears. X av becomes the arithmetic mean in the composed system. 
This is illustrated by the constant for (X\ n ^ al + X^ al )/2. The final average kinetic 

energy per particle, e{*™ 1 = E^"f + {E{ n ^ al + E t 2 n * t v ial )/2, is slightly below the average 
of the initial energies. 

On the figures 3a and 36 we show the equilibration between two systems with different 
non-extensivity parameters: A Boltzmann-Gibbs system (a = 0) and another one with 
a = 2. In 3a the non-extensivity parameter for a collision between particles coming from 
different systems is taken to be a\i = 0. In 36 this parameter is given by a\ 2 = 2. 
In these cases neither X av nor E av approaches to the arithmetic mean of the initial 
X av or E av values (actually X av cannot be defined for the composed system). The final 
common stationary energy distributions are between the initial energy distributions of 
the original subsystems (with ai = and a-i = 2); in case a (a\2 = 0) it is closer to the 
Boltzmann-Gibbs system, in case 6 (a 12 = 2) it is closer to the original non-extensive 
system. 

5. Conclusions 

We studied the equilibration between two non-extensive systems in the framework of a 
relativistic parton cascade model, NEBE. We have demonstrated earlier that these sys- 
tems have non-exponential equilibrium distributions in terms of the free particle kinetic 
energy, which is dependent on the chosen type of non-extensivity (or equivalcntly on 
the details of energy sharing in the two-particle collisions). The composition of different 
non-extensive systems resulted in a common equilibrium distribution. 

We have pointed out that the non-extensive energy addition rule eq. (6), leading to 
a cut power-law distribution in equilibrium, is the next to leading order low-energy 
expansion of a general associative rule. 

We coupled two subsystems after a separate equilibration both by additive and super- 
additive rules. The subsystems tend to achieve a common distribution. This process 
coincides with the equilibration of the temperatures in the case of equal non-extensivity 
parameters. In the case of different non-extensivity parameters, however, neither the 
average of original one-particle energies nor the average of the quasi-particle energies 
is observed. In this case the associativity of the energy composition rule is violated, so 
actually a quasi-particle energy can not be defined. 
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